Magneto-elastic instability in Molecular Antiferromagnetic Rings 
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Lattice stability in a model of antiferromagnetic ring coupled to adiabatic phonons is investigated 
for different values of the spin and numbers of magnetic sites. The magneto-elastic transition is 
shown to be heavily affected by the spin value, displaying a qualitative difference in the nature of 
the instability for spin one-half. Among the different synthesized materials, Cug seems to be the 
best candidate to observe lattice dimerization in these systems. Our analysis excludes stable lattice 
distortions in higher spin rings. The effects of thermal fluctuations are studied in the Cug model, 
where a characteristic crossover temperature is estimated. 
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Magnetic molecular rings (MMR) ,with nearest neigh- 
bors anti-ferromagnetic (AF) interactions, are ideal ma- 
terials for exploring the crossover from microscopic mag- 
netism to collective effects in low dimensional systems. 
Rings with different value of the spin S and number N 
of magnetic sites have been synthesized in recent years 
0. They give the opportunity to understand how the 
static and dynamical properties of the material are af- 
fected by a change in the spin value or a change in the 
number of magnetic sites. Considerable attention has 
been devoted to the observation of resonant tunneling of 
the magnetization in single-molecule superparamagnets 
and of quantum steps of magnetization in molecular 
AF rings H 0. 

clusters, such as 
0, have been studied 
by several techniques, such as nuclear magnetic reso- 
nance (NMR) [3, 13' calorimctric measurements. X-rays 
and inelastic neutron scattering (INS). Recently, INS 
experiments performed on Crg molecular rings ^3 
have been interpreted on the basis of a distorted ring 
geometry at low temperature, in contrast to the cen- 
trosymmetric structure suggested by X-ray diffraction 
at room temperature. In the Cug rings, the Cu NQR 
spectrum shows four structurally non equivalent Cu 
ions An open question is whether the reduced 

symmetry is related to the appearance of dimerization. 
Heat capacity and magnetic torque measurements in 
magnetic field may also suggest the occurrence of a 
dimerized phase in FcgiLi at T = 1 K [lll |. 

The Spin-Peierls (SP) transition is a magneto-elastic 
instability which occurs in a spin chain coupled with the 
lattice. It occurs when a distorted phase characterized by 
a lattice dimerization is stabilized below a critical tem- 
perature Tc due to the gain in magnetic energy. Such 
a transition was first predicted to occur in the (infinite) 
S = 1/2 Heisenberg AF chain and it was indeed 

observed experimentally in the quasi one dimensional 
5=1/2 compound CuCeOg HI. A number of results 
have been obtained both analytically |l4l | and numeri- 



cally ^3 fo'' the spin one-half Heisenberg chain coupled 
to the phonons. Few investigations have be en p erformed 
in the case of finite rings and/or higher spin |lfil |. On gen- 
eral grounds we observe that in the thermodynamic limit, 
the dimer susceptibility of the Heisenberg chain is known 
to diverge as [tt — q\~^ (with logarithmic corrections) for 
wave vectors close to q ~ n and arbitrary semi-integer 
spin S leading to SP instability at T = 0. The 
situation is different in a finite-size system or in chains 
with integer spin 5", where the singlet gap in the exci- 
tation spectrum gives a finite susceptibility also in the 
T ^ limit. Therefore we expect that in MMR the lat- 
tice distortions will be stabilized only for a sufhcicntly 
small elastic constant. The competing role of finite size 
effect and of thermal fluctuations in these systems has 
not been clarified yet. 

In this Letter we first analyze the stability of a meso- 
scopic magnetic ring with respect to lattice distortions 
and then we examine the finite temperature effects in 
these systems. Using Lanczos Diagonalization Tech- 
niques and Spin Wave Theory (SWT) we show how the 
magneto-elastic transition is affected by a change of the 
spin value. Surprisingly, a qualitative difference between 
the behavior of S" = 1/2 and higher spin rings is found, 
strongly suggesting the absence of stable lattice distor- 
tion in the latter case. Among the synthesized MMR, 
Cus seems to be the only candidate to display a magneto- 
elastic instability for realistic values of the phonon cou- 
pling. An estimate of the crossover temperature associ- 
ated to this 'transition" is also given. 

We focus on rings with an even number of sites, in 
which the AF superexchange coupling leads to a singlet 
ground state. We investigate the simplest model Hamil- 
tonian containing both AF interaction and spin-lattice 
coupling 
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where K is the spring constant. Si is the distortion of 
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the bond between site i and i + 1, a is the spin-lattice 
coupling and the last term is the elastic energy. Peri- 
odic boundary conditions are understood (■S'jv+i = ^i). 
Lattice distortions are treated in the adiabatic approxi- 
mation, i.e. Si are c-numbers. Due to the uncertainties 
in the modeling of the microscopic mechanisms leading 
to the spin-phonon coupling in these materials, harmonic 
approximation has been employed in the form of Hamil- 
tonian QJ, thereby limiting its validity to the a\6i\ <C 1 
case. We did not include in the Hamiltonian other terms, 
usually present in the microscopic description of MMR 
[lo| . which do not directly couple with the lattice, such 
as dipolar interactions or spin anisotropics. The invari- 
ancc of the Hamiltonian under the rescaling a Aa, 
K X^K and Si Si/X allows to restrict the parame- 
ter space by setting a = 1 without loss of generality. The 
further constraint, due to periodic boundary conditions, 
J2i = is also understood. By relaxing this constraint, 
the numerical results remain qualitatively unchanged. 

The ground state of Hamiltonian |^ is obtained in two 
steps jlSj : first, by Lanczos technique, we determine the 
lowest energy eigenvalue of at fixed {Si}. The bond 
distortions are then updated in order to comply with the 
equilibrium condition obtained by use of the HcUman- 
Fcynman theorem: 
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J{S, ■ S,+^)s - j;^Y.^S,S,+^)s = (2) 



where (...}<s is the ground state expectation value, in 
presence of distortion. No lattice symmetry has been 
assumed, in order to investigate all the possible lattice 
equilibrium configurations. In all cases we examined, we 
found either an undistortcd configuration or a dimerized 
ring, i.e. Si — (— l)*(5o. No other periodicity has been 
observed, as already known for 5 = 1/2 

The numerical results concerning the stability of the 
undistorted phase are reported in Fig. ^ The distortion 
amplitude is plotted as a function of the elastic constant, 
for different values of the spin. All models we examined 
show a critical constant Kc beyond which the symmetric 
phase is unstable. The results clearly show that, while 
in the spin one-half case the transition is continuous, for 
larger spin the distortion amplitude Sq suddenly jumps 
from zero to a finite value (even larger than the lattice 
spacing) when the spring constant is reduced. Although 
such a large distortion falls outside the range of validity 
of the model, the results clearly point toward a first order 
transition in the S > 1/2 case. No qualitative differences 
between integer and semi-integer spins can be observed 
in these small rings. 

In order to understand the origin of the different be- 
havior of the S = 1/2 case, we have performed Lanczos 
diagonalizations fixing the dimerization amplitude Sq and 
we computed the ground state energy of the dimerized 
Heisenberg model Eh{So)- In terms of this quantity the 
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FIG. 1: Stability diagram for different sizes N and spin S, 
corresponding to synthesized materials. Spring constant K in 
units of J. Left panel: iV = 6 S = 5/2 (Fee); iV = 8 S = 3/2 
(Org) and iV = 8 S = 1. Right panel: = 8, 5" = 1/2 (Cug). 
Inset: blow-up of the Cug stability diagram for K ~ Kc 



equilibrium condition reads 
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(5o = is always a solution, Eh{So) being analytic in Sq. 
Generally, lattice distortions lower the magnetic energy 
Eh for any S, making the undistorted phase stable only 
for K > Kr, with 
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where the prime denotes derivative with respect to the 
dimerization amplitude. The second derivative of the 
ground state energy (i.e. the dimer susceptibility) at 
r.h.s. of Eq. Q can be evaluated numerically for several 
N and S. The resulting critical coupling Kc are shown 
in Fig. 121 suggesting an approximate linear scaling with 
S. Note that the natural energy scale of the problem is 
set by the classical magnetic term ~ JS^. Therefore, 
the dimensionless critical spring constant K^/JS^ ac- 
tually scales as l/S" and vanishes in the classical limit 
S oo. Following the standard Landau-Ginzburg ap- 
proach, the equilibrium distortion amplitude is predicted 
to continuously evolve from zero to finite values as K 
is reduced below Kc. This expectation rests on the as- 
sumption of a positive fourth order term in the expansion 
of Eh {So) about Sq = 0. We checked this condition by 
computing the fourth derivative of Eh at Sq ~ which 
turns out to be positive only in the S = 1/2 case, as 
shown in Table Therefore for S > 1/2, the tran- 
sition is indeed first order as suggested by the numer- 
ical results of Fig. ^ This numerical finding is con- 
firmed by studying the S" — > oo limit of our model 



FIG. 2: Scaling relation for Kc in units of J as a function of 
the spin S, for = 4 (top panel) and N = 6 (bottom panel). 
Black filled symbols are Lanczos data, empty symbols spin- 
wave results 



For a dimcrizcd ring, the classical ground state energy 
Eh = -JUii^ + (-l)"'^o)52 = -^NJS^ becomes inde- 
pendent of the distortion amplitude Jq- It follows that, 
classically, the distortion does not decrease the magnetic 
energy contribution and hence a dimerized geometry is 
not favored. Only quantum fluctuations lift the degen- 
eracy in 5q. In order to estimate the effects of quantum 
fluctuations on the ground state energy we made use of 
lowest order spin-wave theory (SWT) for flnite size chains 
|20l | which is free from singularities even in ID. By ap- 
plying SWT to the dimerized Heisenberg ring, we obtain 
the 0{S) expression for Eh{So)- 



Eh{So) = -NJS{S + 1) + 25 J.\/l - cot (5) 



leading to the result: 
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This explains the linear scaling observed in the numeri- 
cal results. Direct evaluation of the fourth derivative of 
Eh {So) from Eq. (O confirms the first order transition 
in the classical limit 5 — > oo as shown in Table H The 
continuous transition for spin one-half rings is therefore 
exclusively due to the enhancement of quantum fluctu- 
ations and is then a genuine quantum effect which goes 
beyond semi-classical treatments. 

Some further physical insight on the difference between 
spin one- half and S > 1/2 can be gained by investi- 
gating the "proximity" between the undistorted ground 
state and the valence bond (VB) state, corresponding to 
(5o = 1, which mimics the dimerized phase in the strong 
distortion limit. By calculating the overlap between these 
two states shown in Fig. |Hlwe appreciate a noticeable de- 
pendence on the spin value, the overlap getting smaller 
and smaller for increasing S and being significant only 
for S = 1/2. Such a decrease in the overlap suggests that 



FIG. 3: Overlap between the undistorted Heisenberg ground 
state \H > and the VB state \VB > for different spins as a 
function of N 



a continuous transition between the undistorted and the 
distorted phase becomes less and less favored for spin 
larger than one-half at any size N. 

Finally we comment on the stability of the distorted 
phase against thermal fluctuations. The coupling of 
the spins with the lattice (considered here as embed- 
ded in the tree dimensional space) allows for a magneto- 
elastic transition at flnite-temperature, driven by the 
one-dimensional magnetic fluctuations. We concentrate 
on the the behavior of 5 = 1/2 rings for only in these 
cases the ground state may show stable lattice distor- 
tions. In finite size systems no sharp transitions may be 
observed at finite temperature and only crossover tem- 
peratures can be defined. However, a sharp change at 
some Tc{K) is expected in the distortion distributions 
for K < Kc'. The low temperature bimodal distribu- 
tion peaked at ±(5o turns into a Gaussian-like probability 
centered around So ^ for temperatures larger than the 
"mean field" critical temperature Tc{K). By employ- 
ing exact diagonalization in the full Hilbert space of the 
Hamiltonian we computed the partition function and 
the finite temperature dimer susceptibility of the model 
thereby obtaining the critical temperature shown in Fig. 
0] Note the steep increase in Tc{K) close to the transi- 
tion, implying that temperatures higher than ~ 0.1 J are 
always required to destabilize the lattice distortion. 

Now we are ready to apply our results to real mate- 
rials and discuss the possible occurrence of a dimerized 





S = 1/2 


5= 1 


S = 3/2 


S = 2 


S = 5/2 


S = oo 


iV = 4 


54.0 
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TV = 8 


1589.8 


-152.2 


-5.5 






-14.5 



TABLE I: Fourth derivative of the ground state energy with 
respect to the dimerization parameter (5o for different sizes TV 
and spins 5*. Values are divided by 5*. The S = oo results are 
obtained by use of spin wave theory, Eq. 0. 
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FIG. 4: Crossover temperature as a function of the elastic 
coupling for a AT = 8, 5 = 1/2 MMR 

phase in magnetic molecular rings The results of Fig. 
strongly suggest that lattice distortions cannot be ob- 
served in synthesized magnetic rings with S > 1/2. In 
these cases small distortions are never stabilized even at 
T = 0, while the occurrence of a strongly first order tran- 
sition would probably disrupt the ring structure in the 
absence of a sufficiently hard elastic constant K. On the 
contrary, as show in Fie:s. Hll4|l . in the case of Cug {N = 8 
magnetic sites and S = 1/2) a stable dimerized phase is 
instead possible at low temperature. According to Ref. 
0, the Cug ring has a couphng of about J = 0.1 eV, 
the exchange interaction being due to Cu-O-Cu bridge 
as usual in copper-oxide compounds. Assuming a typical 
value for a ~ 10, we find that small lattice distortions 
can be stabilized by a spring constant attaining the phys- 
ically reasonable value K 2± O.lJa^ ~ 1 eV. Instead, in 
iron and chromium rings the magnetic coupling is much 
weaker (J ~ 1.5 meV in Crs ^3 and J ~ 2.4 meV in 
Fee Q 0) pushing lattice instabilities to considerably 
smaller value of K. 

In conclusion, using SWT and Lanczos diagonaliza- 
tions, we have investigated the properties of a generic AF 
magnetic ring, with spins adiabatically coupled to lattice 
distortions. We have emphasized the role of the spin in 
the occurrence of magneto-elastic instabilities, showing 
that the order of the transition strongly depends on the 
value of S. In order to assess the accuracy of the adia- 
batic approximation employed here, experimental studies 
on the characteristic frequencies of the phonon modes will 
be necessary. However, we believe that this test is not 
crucial for the present investigation because lattice insta- 
bilities are classical phenomena whose occurrence can be 



reliably estimated also in the adiabatic limit. Our analy- 
sis suggests that magneto-elastic coupling cannot be in- 
voked for the interpretation of NMR and NQR data in 
MMR with spin different from one-half ^^iJ^^hile con- 
firming this possibility in the case of Cu8 H m. An 
investigation of the effects of thermal fluctuations on the 
stability of the distorted phase has also been performed 
for the latter case. X ray scattering at low tempera- 
ture and NMR spectra analysis will be able to determine 
whether MMR do indeed display lattice distortions in the 
ground state. 
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